Pre-Main Sequence Evolution 


Star formation begins in giant molecular clouds, with temperatures 
less than ~ 100 K or so. In this environment, the average density 
of interstellar matter is ~ 100 cm~, and the size of the cloud is 
anywhere between 10 and 100 pc. Dust is mixed in with the gas, 
which shields the inside of the cloud from UV radiation, keeping 
the temperature cold. (In fact, by measuring molecular absorption 
in these clouds, and assuming a Boltzmann distribution, one can 
obtain the temperature of the microwave background.) 
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The classic picture of star-formation involves a) the formation of 
cores in a giant molecular cloud, b) the creation of proto-stellar 
cores with surrounding disks, c) the breakout of stellar winds, which 
blow material away, and d) the termination of infall, leaving behind 
a star and the remnants of a disk. 




















Figure 7 The four stages of star formation. (a) Cores form within molecular clouds as 
magnetic and turbulent support is lost through ambipolar diffusion. (6) A protostar with a 
surrounding nebular disk forms at the center of a cloud core collapsing from inside-out. 
(c) A stellar wind breaks out along the rotational axis of the system, creating a bipolar flow. 
(d) The infall terminates, revealing a newly formed star with a circumstellar disk. 


The Jeans Criteria 


Before considering the collapse of a proto-stellar cloud, first recall 
the two principal results from our Virial theorem analysis: that for 
a stable star, thermal energy is balanced by gravitational energy, by 


3(y — 1) E; + Egrav = 0 (10.1.8) 
and that the total energy of a star is given by the thermal plus 
gravitational energy. In virial equilibrium, this is 
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where, for a simple monotonic ideal gas, y = 5/3, and (10.1.8) 
reduces to 2E; + Egray = 0. Note that if y < 4/3, the total energy 
is negative, and the system must collapse. 


W=- Poa (10.1.9) 


Now consider the collapse of a gas cloud. Simply put, a cloud will 
collapse if its gravitational potential is stronger than its thermal 
support. (Magnetic and rotational support can also be important, 
but we will not consider those factors here.) A cloud in virial equilib- 
rium with y = 5/3 will have Egrav = —2E;, and, from our discussion 
of the virial theorem, 
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If the gas cloud is uniform, then 
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If we substitute density for radius, i.e., p = M /(4/3)r R°, then 
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This is the Jeans mass, i.e., the minimum mass necessary for a 
collapse. (This derivation is not perfectly correct, since it assumes 
that the gas immediately outside the cloud is static, which is almost 
certainly not the case. It also neglects magnetic support, which can 
be significant.) Note that for a low-density cloud that is optically 
thin, energy is easily radiated. In this very early phase, the collapse 
is roughly isothermal, which means as the density becomes larger, 
the Jean’s mass will become smaller. This leads to fragmentation 
of the cloud and the production of many separate cores. 


Note that if re-write (29.1.4) using the sound speed, cs = (yP/p)!/? 
and assume an ideal, monotonic gas with (y = 5/3), the Jean’s mass 
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Also, we can replace mass with density, to derive the Jean’s length, 
i.e., the maximum size cloud that can collapse 


ik Ne Se Ne 
Ly) ec pee 29.1.6 
sd (a) a (= ( ) 


Another way of looking at the collapse problem is to compare the 
timescale for free fall collapse to the timescale for a pressure way to 
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propagate across the cloud and restore equilibrium. In the absense 
of magnetic pressure, the free-fall collapse time is simply 
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Meanwhile, the time it takes a pressure wave to cross the cloud is 
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Initial Collapse 


As a gas cloud collapses, its increasing density causes it to become 
opaque in the infrared. The trapped energy heats the cloud, and 
changes the nature of the collapse from an “isothermal” phase to 
an “adiabatic phase.” By the definition of an adiabat, 


Pep = T cpl (29.2.1) 
Substituting this into the Jeans equation (29.1.4), means that 
Mj «x p20) pote x gema (20.2.2) 


For y = 5/3, MIX o'/?- in other words, the Jeans mass no longer 
decreases with increasing density. Thus, the cloud is no longer un- 
stable to fragmentation will stop. 

At the center of the cloud, matter is accreted onto the protostellar 
core. As with other accretion, half the energy heats the star, and 
half is radiated away. 
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Since the cloud is relatively faint, the thermal timescale for the cloud 
is large, and the accretion timescale 
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So let’s compute the average luminosity of a protostellar cloud dur- 
ing this phase. In general, for a spherical cloud, the gravitational 
potential energy is 
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Eerav = (29.2.3) 


where n is the polytropic index. Recall that for a uniform density 
sphere, n = 0, and the coefficient becomes the familiar value of 3/5. 
(If one wants to model the cloud as a polytrope with n = 3/2, then 
the coefficient 6/7.) If half the energy is produced during collapse 
of a uniform cloud, then 





1 36M} 
2 5 Ry 





AE ray = (29.2.4) 





Now let’s assume this energy is emitted over one free-fall time, and 
is emitted as blackbody radition (with efficiency, e€). In that case 
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If we replace the Jean’s radius with density, then the minimum mass 
for fragmentation is 


T1/4 
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So, for a typical temperature of T ~ 1000 K, the minimum frag- 
mentation mass is of the order of a solar mass. This is typical of 
star formation: although the initial gas cloud may be hundreds or 
thousands of solar masses, the objects that form are ~ 1 Mo. The 
process is very inefficient! 





Figure 2. The global evolution of the cloud during the calculation. Shocks lead to dissipation of the turbulent energy that initially 
supports the cloud, allowing parts of the cloud to collapse. Star formation begins at t= 1.04tg in a collapsing dense core. By the end of 
the calculation, two more dense cores have begun forming stars (lower left of the last panel) and many of the stars and brown dwarfs 
have been ejected from the cloud through dynamical interactions. Each panel is 0.4 pe (82400 AU) across. Time is given in units of the 
initial freefall time of 1.90 x 10° years. The panels show the logarithm of column density, N, through the cloud, with the scale covering 
—1.5 <log N < 0 fort < 1.0 and —1.7 < log N < 1.5 for t > 1.0 with N measured in g cm~? 





Eventually, the core reaches a temperature of ~ 2,000 K and begins 
to dissociate the molecular hydrogen. At this time, the specific heats 


become very large (since the energy of accretion does not raise the 
temperature of the gas, but just causes a change in the gas state). As 
a result, y = c,/cy < 4/3, and, through the virial theorem (10.1.9), 
the star collapses quickly (on a dynamical timescale). The collapse 
releases energy, which further dissociates H2. Almost immediately, 
all the molecular hydrogen is destroyed, and the star settles into a 
new hydrostatic equilibrium. The same thing happens a short time 
later when the temperature rises enough to ionize hydrogen (and 
then helium, and then helium again). 
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The specific heats and y for a partially ionized plasma with a tem- 
perature of about 10,000 K. In a partially ionized plasma, energy 
goes into increasing the ionization, rather than raising the temper- 
ature. 


Note that a substantial amount of the star’s collapse in this early 
phase is due to these phase-change events. If the proto-star begins 
with a radius Rinit ~ co, then energy associated with dissociation 
(and ionization) is equivalent to 
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where xy, = 4.48 eV, x = 13.6 eV, XHe = 24.6 eV, and XHe+ = 
54.4 eV. In other words, the gravitational energy lost is equivalent 
to ~ 16.4 eV per baryon, and the phase change events alone are 
responsible for bringing a protostar from radius infinity to radius 
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Note also that the temperature required to doubly ionize helium 
(~ 80,000 K) is far less than that required to fuse hydrogen. In 
fact, throughout most of this time, the temperature is cool enough 
so that H~ is the dominant source of opacity. This means that the 
star is fully convective and on the Hayashi track. 





Given that the star is on the Hayashi track (which is essentially 
verticle in the HR diagram), the temperature of this early-evolution 
proto-star is roughly 4,000 K. A size of R ~ 35Rọ then implies a 
luminosity of ~ 10°Lo. Note however, that observed pre-main se- 
quence stars have radii substantially smaller than 35Rọ, indicating 
that additional energy must have been radiated (or otherwise lost) 
in the contraction process. 


Hayashi Track Evolution 


Fully convective stars are located on the Hayashi track. After a 
protostar’s dynamical collapse, it will settle on the Hayashi track at 
a luminosity appropriate to its mass, with an radius approximately 
given by (29.2.8). 


e Evolution will halt on the Hayashi line due to the energy injected 
by deuterium, which burns at ~ 10° K. This energy will up the gas 
pressure and keep the star from collapsing for ~ 10° yr. (Lithium 
burning will also temporarily halt the contraction, but for a much 
shorter period of time.) 


e As the core heats up, more and more of the star becomes dom- 
inated by Kramers law opacity. The star will cease to become en- 
tirely convective (i.e., a radiative zone will develop). When this 
happens, the star will evolve blueward in the HR diagram, towards 
its zero age Main sequence position. 





e As the star approaches the ZAMS, the first steps in CNO pro- 
cessing will occur (even if the core temperature isn’t hot enough 
to support the entire cycle). Specifically, the reactions '*C(p,y)13N 
and 3C(p,y)!4N will convert virtually all the core’s carbon into ni- 
trogen. (This occurs even in low mass stars which burn through the 
proton-proton chain.) This event causes a wiggle in the evolutionary 
tracks as the star approaches the ZAMS. 





e A protostar’s evolution towards the main sequence is on a Kelvin- 
Helmholtz timescale, i.e., 
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Since £ is a high power of mass, the time it takes a star to reach 
the main sequence depends strongly on mass. Moreover, as the 
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star contracts, R gets smaller, which implies Tx 7 gets larger. The 
star’s pre-main sequence lifetime is therefore dominated by the final 
stages of contraction, when the star is already close to the ZAMS. 
For stars with M > 1 Mọ, TKH ~ 5 x 107(M/Mo) 7? yr. 
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Note that the pre-main sequence timescale implies that in young 
star clusters, high mass stars will evolve off the main sequence before 
low-mass stars have finished contracting! 
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Identification of Young Stars 


There are several ways to identify Young Stellar Objects: 





e By infrared excess: Young stars are built up by accretion, and the 
accretion disk can linger for some time after the nuclear reactions 
have started. The dust and material in the disk can reprocess some 
of the light into the IR. 
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e By X-ray identification: Conservation of angular momentum guar- 
antees that young stars will be rotating quickly. If the star has 
a convective envelope, then differential rotation will cause an in- 
creased magnetic dynamo effect, leading to increased flare activity 
and X-rays. 
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e Ha emission: The same mechanism that can create X-rays may 
also result in Ha emission from prominences, etc. In addition, there 
may still be residual Ha emission from the accretion disk. 


e Lithium absorption: Lithium can be burned during the pre-main 
sequence phase. Moreover, in stars with convective envelopes, the 
temperature at the bottom of the convective layers is sufficient for 
lithium burning. Stars with lithium absorption in their spectrum 
must be young. 
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